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A Higgs boson produced in association with one or more bottom quarks is of great theoretical and 
experimental interest to the high-energy community. A precise prediction of its total and differential 
cross-section can have a great impact on the discovery of a Higgs boson with large bottom-quark 
Yukawa coupling, like the scalar (/i" and H^) and pseudoscalar (A") Higgs bosons of the Minimal 
Supersymmetric Standard Model (MSSM) in the region of large tan/3. In this paper we apply the 
threshold resummation formalism to determine both differential and total cross-sections for bg — > 6$ 
(where $ = h^,H"), including up to next-to-next-to-next-to- leading order (NNNLO) soft plus virtual 
QCD corrections at next-to-leading logarithmic (NLL) accuracy. We present results for both the 
Fermilab Tevatron and the CERN Large Hadron Collider (LHC). 

PACS numbers: 13.85.-t, 14.80.Bn, 14.80.Cp 



I. INTRODUCTION 



^ ' The prospect of discovering the mechanism of Electroweak Symmetry Breaking (EWSB) in the coming years is of 
, primary interest to high-energy physics programs around the world. In the Standard Model (SM), a single SU{2)l 
complex scalar doublet gives rise to one Higgs boson (h) and massive gauge bosons through the Higgs mechanism as 
well as massive fermions through Yukawa interactions. In the Minimal Supersymmetric Standard Model (MSSM), 
^ . there are two complex scalar doublets each giving rise to the masses of the up- and down-type quarks respectively. 

' After EWSB in the MSSM, five physical Higgs bosons remain: a light and heavy scalar (/i°, H^), a pseudoscalar 
CO ■ and two charged scalars {H^). Their couplings to the SM particles may substantially differ from the SM Higgs boson 

■ couplings. Many of the properties of the SM and MSSM Higgs bosons have been reviewed in Refs. [H, 

I There have been some rapid changes in the precision electroweak fits that have lead to new bounds on the Higgs 

■ boson mass in the SM, mostly from a shift in the central value of the top quark mass 3. Together with the 
CD ] lower bound set by direct searches for a SM Higgs bosons at LEP-2 precision electroweak fits indicate that 

'. 114.4 < Mh < 166 - 199 GeV at the 95% confidence level The experimental bounds on the MSSM Higgs bosons 
are weaker than those of the SM due to the much larger parameter space of the supersymmetric model, leading to 
a Mf^o j^o > 93 GeV lower bound 0, However, the lightest MSSM scalar is theoretically bounded to lie below 
about 130 GeV. Therefore at least one Higgs boson that may exist in nature (either SM or MSSM) will probably be 
^ ,, constrained by the Tevatron and definitely probed by the Large Hadron Collider (LHC). 

In the SM, for current and future planned energies at hadron colliders (both pp and pp), a Higgs boson would be 
primarily produced via gluon fusion {gg — > h) through a top-quark loop. However, in the MSSM the production 
channels can become quite varied depending on the parameter space. In particular, for large values of tan/3 (the ratio 
of the two Higgs doublets vacuum expectation values) the MSSM bottom-quark Yukawa couplings are enhanced and 
the production of a Higgs boson with bottom quarks become the leading production mode. Indeed, this has already 
been used at the Tevatron to substantially constrain the MSSM parameter space Q. 

Bottom quarks pro duced in association with a Higgs boson have also raised a great deal of theoretical interest 
Hoi , [ill . [T^ . [isl . [l^. llSl. [l6l . [iTj . Although there is a conceptual difference in whether one considers bottom quarks in 
the initial state as partons or only gluons and light quarks (i.e., whether one employs the five fiavor number scheme, 
5FNS, or the four flavor number scheme, 4FNS), physical observables have been shown to agree remarkably well in 
both schemes once full next-to-leading order (NLO) QCD corrections are included [H, [13, [lal- This has placed the 
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theoretical predictions of both total and differential cross-sections under better control. Further improvement can be 
achieved by considering the impact of resumming corrections that can be relevant in specific regions of phase space. 

In this paper we investigate the impact of higher-order QCD threshold corrections on the bg &i> ($ = total 
and differential cross-sections, both at the Tevatron and at the LHC. The resummation of soft plus virtual dynamical 
threshold corrections to a variety of p rocesses involving both scalars and pseudoscalars Higg s bosons has been studied 
extensively in the literature [l3, [3 H [lU, M, M, M, MM^^MiM, El, M, M, M, iHf . In particular, following 
the resummation techniques originally proposed in Refs. 34. l35l~36l|. we will use the formalism recently outlined in 



Refs. |32, i38i, |39| , to which we refer for further references. 

In order to assess the validity of the threshold resummation formalism for the case of bg = h^,H^) 

production, we compare the NLO full calculation (in the 5FNS) Q as obtained from MCFM [l^l with the results 
obtained by expanding the resummed cross-section to include NLO soft plus virtual corrections up to next-to-leading 
logarithmic (NLL) accuracy. The NLO-NLL formalism reproduces very closely the full NLO calculation in the region 
of small Higgs boson transverse momentum, where most of the statistics are accumulated. We then improve upon the 
NLO fixed-order calculation by including both next-to-next-to-leading (NNLO) and next-to-next-to-next-to-leading 
order (NNNLO) soft plus virtual corrections up to NLL accuracy. Higher-order corrections are more and more 
stable with respect to variations of both the renormalization {(ir) and factorization (fip) scales, in particular in the 
threshold region where soft corrections dominate the bulk of radiative corrections. We will discuss more extensively 
the perturbative behavior of NNLO-NLL and NNNLO-NLL corrections in Section ITVl 

Overall, the NNLO-NLL and NNNLO-NLL effects are sizable and they greatly stabilize the cross-section. The 
improved small transverse momentum behavior should help to shed some light on the current bounds on tan /3 in a 
more robust way, since it corrects the distribution in a region where most of the statistics are accumulated. 



II. MSSM PARAMETERS 



Higgs boson physics when associated with bottom quarks is dominated by two factors - the running of the bottom- 
quark mass and the MSSM couplings to the different CP-even and CP-odd neutral Higgs bosons. The one-loop 
supersymmetric corrections to the bottom-quark Yukawa coupling must also be considered. We will briefly consider 
each of these factors to establish our notation. 

The SM bottom-quark coupling to the Higgs boson is g^^^ = rnb{fJ.R)/v, where the vacuum expectation value (v) 

of the Higgs doublet is defined as v = (^/2G^^)-l/2 = 246 GeV and mb(^i^) is the MS running mass of the bottom 
quark as a function of the renormalization scale ^r. We define rnt,{fiR) up to 4 loops (as we will need for the NNNLO 
corrections) according to Ref. [4l[, which corresponds at 1 and 2 loops to the following expressions: 
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where a^diR) is the strong coupling constant at the renormalization scale fiR and. 
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are the first two coefficients of the QCD beta function and mass anomalous dimension function, with i\f = 3, the 
number of colors, and = 5, the number of light fiavor. The quantity is the bottom-quark pole mass which 

we take to be m^°''' = 4.62 GeV. The 3-loop and 4- loop expressions of mi,{fiR) are lengthy and we refer to Ref. [4lj 
for their exact definition. In the following, we rescale the 5FNS NLO results obtained from existing calculations via 
MCFM in such a way to agree with this definition of rni,{fiR). This is possible because Tni,{fj,R) enters only in the 
overall bottom-quark Yukawa coupling. 

The tree level CP-even neutral Higgs boson couplings to bottom quarks in the MSSM can be written in terms of 
the SM couplings as. 
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and are enhanced in the limit of large tan /3, where tan(/3) = V2/V1 is the ratio of the vacuum expectation values of the 
Higgs doublets coupling to the up- and down-quarks respectively, while a is the angle which diagonalizcs the neutral 
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Higgs sector of the MSSM. The CP-odd neutral Higgs boson couphng to bottom quarks is also enhanced exactly by 
a factor of tan/?. In the following, we will focus our discussion on the CP-even neutral Higgs boson cases. In the 
mi, ~ limit, the results obtained for the neutral scalars can be directly rescaled to obtain results for the neutral 
pseudoscalar. For non-zero mf,, the difference between the scalar and pseudoscalar case, modulus a rescaling of the 
couplings, is of order (mf,/Af$)^ {^ — h'^,H'^), and one therefore expects it to be small. 

Finally, we need to consider the supersymmetric corrections to the bottom-quark Yukawa coupling in Eq. ^ due 
to squark and gluino loops. For CP-even neutral Higgs bosons we find [4^ . 
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where the one-loop correction can be written as, 

2as{mt) 
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and the integral quantity /(a, 6, c) is defined as. 



I[a,b,c) 



(a2-62)(62_c2)(Q2_c2) 



(4) 
(5) 

(6) 

(7) 



Other MSSM parameters of influence in the previous set of equations are the Higgs-Higgs coupling in the super- 
potential, /Li, the masses of the up- and down-type squarks after mixing, m^^ ^ and m^^ ^, the mass of the gluino, nig, 
the SM top-quark Yukawa coupling gt^j^ = rrit/v, and the top-quark tri-linear coupling, Af . We set nit = 172.2 GeV, 
mg^ ^ = ^ = MsusY = 1 TeV, rUg = 1 TeV, A^ = At = 2 TeV, and fi = M2 = 200 GeV, where M2 is the 
SU{2) gaugino mass parameter. With these values we find Af, = 0.178. In our convention, the sign of Higgs-Higgs 
coupling, n, is preferred positive. All of the MSSM couplings were calculated using the FeynHiggs 2.5 package [43| 
which includes all known corrections to the bottom-quark couplings in the MSSM through 2 loops. 



III. RESUMMATION 



There are several formalisms for resummation calculations. Much of the research in resummation focuses on total 
cross-sections, however there are several ways of approaching differential cross-sections as well [2TI. [28l. [32l. [33| . 

Different formalisms lend themselves better to different processes. In general terms, for differential quantities, we 
can consider the resummation of one particle inclusive (IPI) or pair invariant mass (PIM) observables. In the case of 
our observable, we are interested in resumming a 2 ^ 2 process, bg — s- 6$ {^ — h'^^H'^), so we are in the domain of 
the IPI dynamics. The formalism that we will use to resum threshold corrections to the differential cross-section of 
a 2 — > 2 process, implementing known universal corrections, is due to Kidonakis [3^. Issl. Is^ . 

The calculation of both total and differential cross-sections in hadron-hadron collisions can be formalized as 

0- = ^ dxi dx2(|)i/h^{Xl,^iF,^^R)4'j/h2ix2,fJ.F,fJ■R)^iJis,ik,fJ.F,fJ■B.) , (8) 

t.J 

where a and Uij could be a total or differential cross-section of interest, at the hadron and parton level respectively, 
while the indices i and j run over the species of partons contributing to a given process, ipi/hi (for I = 1,2) is 
the Parton Distribution Function (PDF) for parton i carrying a fraction xi of the momentum of hadron /i;, at a 
factorization scale jjLp and renormalization scale /i^. The partonic center of mass energy is s and tk are partonic 
t-channel type Mandelstam variables. 

Higher-order QCD corrections to total and differential parton-level cross-sections often contain so-called plus dis- 
tributions and delta functions of the kinematic variables that in essence measure the deviation from the kinematic 
threshold. If we define S2 in terms of the partonic Mandelstam variables SiS.S2 = s + t-\-u — mf, which vanishes 
at threshold, we will find plus distributions 2?; (§2) and delta functions 5{s2) in the differential cross-section where, 

S2 _ , 
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and M denotes a typical mass-scale of the process. At order a", we expect Z < 2n — 1 plus distributions. When 
the differential or total cross-sections are resummed in moment space, these plus distributions appear as logarithmic 
divergences in terms of the moment variable. These are the logarithms that are resummed at all orders by the 
threshold correction resummation formalism. Close to threshold, i.e. in the region where any extra emitted parton is 
necessarily soft, these corrections can be large and even dominate the cross-section. Hence the need to resum them. 

The formalism in Ref. [13, HH allows us to write these resummed corrections in terms of the tree level cross-section 
times universal coefficients based on the color ffow of the process. The color flow of our process, hg &$ = 
is the same as charged Higgs production, hg H^t, that has been studied in Refs. ^37, jig, 39j. 

The resummed cross-section cannot be easily evaluated outside of its moment space [331 ' however, it can be expanded 
and evaluated in the usual manner. Expanding the resummed cross-section to include NLO, NNLO and NNNLO soft 
QCD corrections at the NLL accuracy provides us with the expressions we will use in this paper. According to 
Refs. [13, [H, [1^, the corrections specific to our process can be written in terms of three coefficients, ci,C2, and C3, 
determined by the color ffow and the kinematic invariants of t he p rocess. For completeness, let us explicitly give ci, 
C2 and C3 in the following equation, using the notation of Ref. [38|: 
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where T, 



2 , the scale-independent part of the c'^^^'^ coefficient, is defined by 
rj^bg^b^ = 2Rer'j'' ~Cf~Ca- 2Cf In 
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and r'^^', related to the one-loop soft anomalous dimension which describes non-coUinear soft gluon emission, is given 
by: 
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In Eqs. (Unil-dni) we have used Ca^Nc^Z, Cf^{NI - l)/27Vc = 4/3, /3o = (HCa - 2n/)/3, n/==5, while M^ is the 
mass of the neutral CP-even Higgs boson (<!> = i/°) and s, t and u are the parton level Mandelstam variables for 
the process b{pi)g{p2) &(— P3)3'(— Ps), defined as s={pi +P2)^, t={pi +^3)^, and u = {p2 +^3)^- The other mass 
scales in this process include p§ = 



fig— >&$ 



;^ and = Adl- We notice that, following Ref. [38|, we have included in c{ 
only the scale-dependent pieces of the S{s2) corrections. In particular, we do not include in cj^~*^* the full virtual 
corrections. We have estimated the impact of the scale-independent terms using the full one-loop virtual corrections 
calculated in Ref. [^ in the mb = limit, and we have found them to be very small. 

In terms of the ci , C2 , and C3 defined in Eq. ITUl) , the NLO-NLL corrected (total or differential) cross-section can 
be written as: 



— {C3I?1(S2) + C2Vo{s2) + Ci6{§2)} 
TT 



(13) 



where o-^™" is the parton level tree level cross-section (total or differential) for bg 
level differential cross-section for bg — > 6$ can be written as: 



6$. For completeness, the tree 



di du 



V ) 

{-t[il4 + u^] + m^[4M| - s2 _ 2u{Ml + s)] - m^[6s + Zt + Au]] 



NciK^ - 1) 
1 

s(t — ?Tl^)2 



(14) 



where the Mandelstam variables in the partonic system were defined above and 



pole 



. We notice that we have 



treated the final state bottom quark as massive in the bg — > &$ tree level cross-section. The result in Eq. ([T?|) matches 
the well known result [^ in the limit of vanishing bottom-quark mass. There are many reasons to keep the final 
state bottom quark massive in this calculation. The prcfactor in Eq. (|14p clearly shows that the bottom quark mass 
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regulates the small transverse momentum region of the born level differential cross-section, allowing integration down 
to zero transverse momentum - an important improvement in our treament of this process. The initial state bottom 
quark is left massless to be consistent with its treatment in the evolution of the parton distribution functions. This is 
consistent with current treatment of the heavy flavor thresholds [3 IH, IH, 13 '^^^ chosen PDF sets and is also 
consistent with the implementation of the splitting functions for the energy evolution which are all strictly massless 
and assume massless initial state quarks. 

In the same way as the above differential cross-section, the NNLO-NLL corrected (total or differential) cross-section 
can be written as [331, 
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while the NNNLO-NLL corrected (total or differential) cross-section can be written as [3£ 
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where C2 = C(2) = 7r^/6, Cs = C(3) = 1.2020569..., and 4 are the scale dependent parts of ci and C2 defined in 

2 > 

Eq. IIOI), r2 is defined in Eq. (HU^hile X^^f , Xf , and X-- (for i = 1, 2, 3) are functions of the kinematic variables 



of the process described in Refs. [38| and [3£ 

It is important to understand the limitations of this calculation. Indeed, the kinematics of the rcsummcd total or 
differential cross-section is fixed by the tree-level process (see Eqs. ([T3|)-([IH)) and it is therefore a 2 ^ 2 kinematic, 
even upon inclusion of resummed higher-order QCD corrections. At the same time, the observation of Higgs boson 
production with one 6-jet is subject to identification cuts imposed on the final state 6-jet transverse momentum (p^) 
and pseudo-rapidity (r^fc), which translate, given the 2 — > 2 kinematic of the resummed approach, into cuts on p*, 
at all orders of the resummed QCD corrections. This has to be taken into account when looking at distributions, 
in particular the p* distribution. Therefore, we interpret the corresponding distribution as an improved estimate 
of the corresponding fixed-order result in the region of low p* above the cut. If the experimental cut on is 
lowered, larger and larger portions of the resummed distributions become important, since one enters more and more 
the region where the resummed corrections are important and a fixed-order calculation is not expected to give a 
reliable result. In spite of the fact that the resummed p* distribution does not include effects from gluon dynamics, 
the bottom-quark mass does much to stabilize the small transverse momentum region and gluon dynamics effects 
should be under control in this region. 

There is an additional complication when we compare our resummed results to the fixed-order results produced by 
MCFM. The way in which MCFM calculates the NLO corrections to Higgs plus bottom quark processes is broken 
into three parts, as described in the MCFM manual [40|. The process labeled 143 takes into account a Higgs produced 
with two additional bottom quarks, both of which are observed. This process would need to be treated separately in 
our resummation formalism. The unaccounted for process has physical significance for our signal if at least one of the 
bottom quarks were to meet the identification cuts. Therefore we have added this process to our resummed results. 
This addition makes our hybrid curves highly reliable in all regions. 
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IV. RESULTS 

In this section we present several results that will illustrate both formal and phcnomcnological aspects of the 
rcsummed calculation. We always show results for bh^ production at the Tevatron and for production at the 
LHC. The final state bottom quark is identified imposing that its transverse momentum and pseudorapidity satisfy: 
> p^*^"* (where p^'^"* will be specified separately for different plots) while \r]''\ < 2 at the Tevatron and jry^j < 2.5 
at the LHC. We will specify when no identification cuts arc imposed on the rcsummed calculation and further justify 
our choice. 

With the exception of Fig. [3l where we study the rcnormalization and factorization scale dependence of the 
rcsummed cross-section, we set the renormalization and factorization scales to fij^ = /ij? = fio/2, where fiQ = 
(ml"^" + M$/2)[5|]. LO results are obtained using CTEQ6L1 PDFs 1-loop as{fJ,R) and 1-loop bottom-quark 
running mass 

(TO6(Mi?)) (see Eq. (P)). NLO results are obtained using the NLO set of MRST 2004 PDFs |50l (since 
they also provide NNLO PDFs), 2-loop asifin) and 2-loop rnb{nR) (sec Eq. (P)). Finally NNLO and NNNLO resuhs 
are obtained using the NNLO set of MRST PDFs [lOl, as well as 3-loop and 4-loop as(^fl) and rnb{fiR) [4l[ respec- 
tively. Although the MRST group has produced a LO PDF set [sij, it is from an older fit to data than the CTEQ6L1 
PDF set. We also note that the available LO MRST 2001 set includes NLO and NNLO PDFs which differ from the 
MRST 2004 NLO and NNLO sets, as expected, and keeping a consistent PDF set is no longer viable option if modern 
PDFs are required. Moreover, as our lowest order results are meant for strictly illustrative purposes, the choice of 
modernity over uniformity is a minor one. 

Figs. [1] and [5] arc obtained using MSSM couplings, since they arc of direct phcnomcnological interest. In particular 
we use the setup explained at the end of Section [TTl with tan/3 = 40. Different M^fi fjo hence correspond to different 
values of Mao. On the other hand, the perturbativc properties of the rcsummed cross-sections (see Figs. [21 [3] and |4]) 
are studied using SM couplings. 

In Fig.[T]we plot the p* (<& = /i", H^) distribution for Higgs production with one bottom quark at both the Tevatron 
and the LHC, for A/^o = 120 GcV and A/^o = 200 GcV respectively. We use MRST 2004 PDFs for both fixed-order 
and rcsummed results. The fixed-order results arc obtained using MCFM [5^. The NLO-NLL and NNLO-NLL 
resummation results, obtained with no p^ identification cut imposed on the final state bottom quark, arc compared 
to the NLO 5FNS fixed-order results with decreasing values of the p^ identification cut (p^ > 10,5,2, and 1 GcV). 
As we already commented at the end of Section Hill the resummed cross-section in Eqs. (|T^ - (|TB1) is bound to a 2 ^ 2 
kinematic. A cut on p\, automatically translates into a cut in p* and truncates the p* spectrum below that point. 
Therefore, in order to obtain the p* distribution over the entire p* range, we do not impose a cut on the py of the 
resummed differential cross-section with the caveat that fixed-order distributions obtained for a given cut on p^ have 
to be compared with the rcsummed distribution for p* > p^'^"*. 

Overall, even for sizable cuts on p^ (say py > 10 GeV), we can see very good agreement between the NLO fixed- 
order calculation and the NLO-NLL resummed differential cross-section well below the mass scale of the process [M^o 
or Mfjo), and we begin to see an appreciable difference starting at ^r = fip = /^o/2 = (m^°'''-|-M$/2)/2, as expected. 
The resummed distribution seems to interpolate between the fixed-order curves, smoothing the low p* portion of 
the spectrum where the large soft contributions start to be more relevant. Higher-order terms in the perturbative 
expansion of the resummed cross-section should then represent better approximations of the p* spectrum in the low 
momentum region, where most of the statistics are accumulated. We then interpret the NLO-NLL and the NNLO- 
NLL curves as improvements over the fixed-order calculation in the region of low p* above the experimental p^ cut. 
For the choice of parameters (Higgs boson masses and renormalization/factorization scale) and PDFs in Fig. [T] the 
NLO-NLL and NNLO-NLL curves are remarkably close. For other choices of parameters they typically have the same 
shape, but they can have more distinguishable magnitudes. This is better illustrated in Fig. [31 where we study the 
dependence on the rcnormalization and factorization scales of the various terms in the expansion of the resummed 
cross-section. 

We have also compared the total cross-sections for bh^ and bH'^ production when a p^ identification cut is imposed 
on both fixed-order and resummed calculations. In the resummed case, this corresponds to dropping the region of 
the p* spectrum below p^. By looking at Fig.jTlwe do not expect the fixed-order and resummed results to agree well 
when large py cuts are imposed. Results are shown in Fig. [2l where the total cross-section is plotted as a function of 
($ = h°, H°), for p^ > 10 GeV and p^ > 40 GeV respectively, at both the Tevatron and the LHC. The Tevatron 
plots, in particular, confirm our expectations. We also notice that the NNLO-NLL and NLO-NLL results converge 
better the larger the Higgs boson mass, i.e. when the kinematic approaches the threshold condition. This last aspect 
of the perturbativc behavior of the resummed cross-section is strongly confirmed in Fig. [3l 

Indeed, Fig. O shows the perturbative behavior of the resummed cross-section expanded at LO, NLO-NLL, NNLO- 
NLL and NNNLO-NLL, when the renormalization and factorization scales are varied between fio/2 and 2/io, for 
jiQ — m-b"''' + A/$/2 ($ = /i°,_ff°). We keep hf = Mi? in our variation since this is how PDF packages arc structured 
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FIG. 1: Comparison of the NLO fixed-order 5FNS and the NLO-NLL and NNLO-NLL resummed transverse momentum 
distribution at the Tevatron and the LHC for different cuts in the transverse momentum of the final state bottom quark in 
the fixed-order calculation. We have added the additional final state 6b<& ($ = h",!!'^) to our resummed results as described 
in the text. The arrows at the top of the graph guide the eye to show where each of the fixed-order distributions peak 
when the histogram is smoothed. As expected, lower cuts on the transverse momentum lead to peaks at smaller values. 
The fixed-order differential cross-section also begins to pick up more of the large soft contributions from the region of small 
transverse momentum. Here we have set M/^o = 120 GeV (Tevatron), M^o = 200 GeV (LHC), and hr = /if ~ /io/2 with 
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FIG. 2: A comparison of the resummed NLO-NLL and NNLO-NLL total cross-sections with the fixed-order NLO corrections 
as calculated by MCFM (5FNS), for > 10 GeV and > 40 GeV, with varied Higgs boson masses. We have added the 
additional final state fe6$ ("I> = h^\H^) to our resummed results as described in the text. 



and varying the two scales separately would cause a mismatch at this level. The LO, NLO-NLL and NNLO-NLL bands 
are completely consistent, since all pcrturbative quantities {as{^lR)^ rnb{iJ-R) and PDFs) are defined at the correct 
perturbative order. The NNNLO-NLL cannot be matched to the correct perturbative order PDFs, since they are not 
available at NNNLO. No identification cut has been used in this plot, since we are only interested in the theoretical 
behavior of the cross-section as a function of the scale. For the same reason, we have used SM couplings. Wc notice a 
large impact in going from LO to the first order of QCD corrections, i.e. NLO-NLL, while the NNLO-NLL order add 
only a small variation. The large difference in going from LO to NLO could in part be due to the different set of PDFs 
used. On the other hand, within the theoretical uncertainty due to the residual scale dependence, the NLO-NLL and 
NNLO-NLL predictions are completely consistent. The residual uncertainty in the NNLO-NLL prediction is smaller 
than the corresponding uncertainty in the NLO-NLL prediction for large Higgs boson masses (see inlays in Fig. ([3])), 
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FIG. 3: A comparison of different perturbative orders of the resummed total cross-sections versus Higgs mass when the 
factorization fip and renormalization fin scales are varied by a factor of two about the central value hf ~ f^R = A'o with 
^0 = mP°''= + AU/2 ($ = h°,H"). Each band shows the maximum variation for a given perturbative order. The inlays 
illustrate in greater detail the large M<i, region. We have added the additional final state 66$ ($ — h",H'^) to our resummed 
results as described in the text. 



i.e. in the threshold region, where the resummed formahsm works better. 

All of these characteristics neatly illustrate the correct behavior of the resummed perturbative series. As far as the 
NNNLO-NLL corrections go, they do not seem to follow this pattern of compatibility with the lower order result and 
reduction of the residual scale dependence, but we need to remember that there is a mismatch with the perturbative 
order of the PDFs. Because of this, we would suggest to limit the improvement of the NLO fixed-order cross-section 
to the NNLO-NLL resummed predictions. This is indeed what we show in the distribution plots of Fig. [TJ 

In Fig. |4l we summarize our results by plotting various K-factors, i.e. ratios of the higher-order (NLO-NLL, 
NNLO-NLL, and NNNLO-NLL) to the LO cross-section, as a function of the Higgs boson mass. We also give ratios 
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FIG. 4: Ratios of the NLO-NLL, NNLO-NLL, and NNNLO-NLL to the LO total cross-section for different values of the Higgs 
boson mass (M^o at the Tevatron and Mjjo at the LHC), and fj.R = /if = ^o/2 = {rn^°^'' + Af/iO,jjo/2)/2. We have added the 
additional final state 66$ ($ — h}\ H°) to our resummed results as described in the text. 

of the various orders of the resummed cross-section to the NLO fixed-order cross-section, calculated with MCFM. AU 
higher-order results use MRST 2004 PDFs, while the LO ones use CTEQ6L1. All curves on Fig. 2] are obtained for 
I^R = /i_F = /io/2 for yUo = + These curves quantify the effect of the higher-order corrections to the total 

cross-section, as well as to the fixed-order cross-section. We see that the threshold effects are more relevant at the 
Tevatron than the LHC for these particular parameter choices, in that the ratios deviates more from unity at the 
Tevatron than at the LHC. 

Finally, in Fig. O we emphasize the behavior of the bh^ and hH^ cross-sections in the MSSM, by plotting both 
production cross-sections in the same plot as a function of M$ ($ = /i", The specific behavior of the MSSM Higgs 
boson masses and couplings is very well represented. As the lightest Higgs scalar boson {h^) becomes inaccessible, the 
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FIG. 5: Total MSSM cross-sections for the two processes, bg ghP and bg — > faff", at the Tevatron and the LHC respectively, as 
a function of Mf^ jjo , when iir = fip = ^J^o /"^ = ('Ti^°'°+i\fftO /2)/2. We have added the additional final state (<1> = h°,H'^) 
to our resummed results as described in the text. 

heavier scalar Higgs boson turns on and begins being produced. The effects are on top of the QCD corrections, 
that are here represented in terms of the various terms in the perturbative expansion of the resummed cross-section. 
The results plotted in this figure have been obtained without imposing a specific identification cut on the final state 
bottom quark, since the plot only aims at showing at a glance the effect of adding leading higher-order terms in the 
calculation of the total cross-section. It is clear that all QCD effects should be included in any realistic attempt to 
find a Higgs boson(s) in nature in the coming years. 

V. CONCLUSIONS 

We apply the threshold resummation formalism to improve the NLO fixed-order predictions for the production 
cross-sections of a MSSM Higgs boson {h'^, H'^) in association with one bottom quark. We focus in particular on the 
Higgs boson transverse momentum distribution and study how to use the the threshold resummation formalism to 
provide a more accurate prediction in the region of low transverse momentum. We also study in detail the perturbative 
behavior of the resummed cross-section and establish its limits and validity. 

We see that the threshold corrections are important at both the Tevatron and the LHC, in particular for large Higgs 
boson masses, because the Higgs would be produced closer to threshold where the resummation effects are greatest. 
As expected, at both colliders the most relevant impact is in the low portion of the distribution, where also most 
of the statistics are accumulated. Therefore, we expect these results to provide valuable information to experiments 
searching for evidence of MSSM scalar Higgs bosons at both the Tevatron and the LHC. 
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